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Abstract 



The diffractive dissociation of 500 GeV/c vr" into di-jets is described as a 

Q . way to measure the momentum distribution of quarks in the pion. The mea- 

^j I surements of the pion diffractive dissociation were carried out using data 

^f^ ' from Fermilab E791. Preliminary results show that the \qq) Asymptotic 

wave function which was developed using perturbative QCD methods de- 
scribes the data well for Q^ ~ 10 (GeV/c)^. At these values signals of color 
^ ' transparency are expected and, indeed, observed through the ^-dependence 

■^ ■ of the yield of the diffractive di-jets. 
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a^ ■ 1 Introduction 

^ I One of the most sensitive tests of QCD is measurement of the internal momentum 

I I distribution of quarks in hadrons. These distributions were calculated for the 

Q^' pion already about 20 years ago (section 2.1) through calculations of the light- 

^ ■ cone wave functions. However, there were no direct experimental measurements. 

Measurements of observables which are related to these distributions, such as 

k> ! the pion electromagnetic form factors, turned out to be rather insensitive to the 

^ I light-cone wave functions. The pion wave function is expanded in terms of Fock 

states: 

^ = a\qq) +(3\qqg)+-f\qqgg) + .... (1) 

The first (valence) component will be dominant at large Q^ as the other terms are 
suppressd by powers of w? /Q"^. 

2 The Pion Momentum Wave Function. 

2.1 Theoretical predictions 

Two functions have been proposed for the \qq) configuration. The Asymptotic 
function was calculated [1, 2] using perturbative QCD methods and is thus ex- 
pected to be correct for very large Q^ (Q^ -^ oo): 
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X is the fraction of the longitudinal momentum of the pion carried by a quark in 
the infinite momentum frame (see Fig. 2). Using QCD sum rules Chernyak and 
Zhitnitsky [3] proposed a function that is expected to be correct for low Q^: 
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The wave functions are plotted in 
Fig. 1 and, as can be seen, there is big 
difference between the two. Measure- 
ments of the electromagnetic form fac- 
tors of the pion were considered as the 
best way to study these wave functions. 
A comprehensive summary of the sta- 
tus of these measurements was recently 
published [4] . Measurements of the pion 
electric form factor suffer from two ma- 
jor drawbacks: measurements done by 
elastically scatterred pions from atomic 
electrons measure directly the pion elec- 
tric form factor but can be done only at 
very low Q^. The alternative method 
is measurements of p{e, eV^)n in the 
electron-pion quasi-free scattering kine- 
matics. The relation between the form factor and the longitudinal cross section 
is model dependent as it includes the p -^ mr^ matrix element. Finally, the form 
factor is related to the wave function through integral over the wave function 
and the scattering matrix element reducing the sensitivity to the wave function. 
Indeed, as shown in [4] the two wave functions can be made to agree with the 
experimental data. Similar situation exists for inelastic form factors, decay modes 
of heavy mesons etc. The problem is that comparisons with these observables are 
not sensitive to details of the wave function and thus cannot provide critical tests 
of the a;-dependence of the functions. An open question is also what can be con- 
sidered as high enough Q^ to qualify for perturbative QCD calculations and what 
is a low enough value to qualify for a treatement based on QCD sum rules. An 
experimental study that provides information about the momentum distribution 
of the \qq) in the pion will be described. 



Figure 1: Two predictions of the wave 
functions of the pion. Asymptotic func- 
tion: dashed line, and CZ function: 
solid line. 



2.2 Pion Diffractive Dissociation and Momentum Wave 
Function. 

The concept of an experiment that maps the pion momentum wave function is 
shown in Fig. 2. A high energy pion dissociates diffractively on a nuclear target 



XP:, 



(i-x)P» q 



t 



imparting no energy to the target so that it does not break up. This is a coherent 
process in which the quark and antiquark break apart and hadronize into two jets. 
If in this fragmentation process the quark momentum is transferred to the jet, 
measurement of the jet momentum gives the quark (and antiquark) momentum. 
Then: 

_ Pjetl / -N 

•^measured , \^) 

Pjetl + Pjet2 

By varying kt (Fig. 2) the measurement can map the x dependence for different 
\qq) sizes. 

The diffractive dissociation of high 
momentum pions into two jets can be 
described, hke the inclusive Deep Inelas- 
tic Scattering (DIS) and exclusive vector 
meson production in DIS, by factoring 
the perturbative high momentum trans- 
fer process from the soft nonperturbative 
part [5]. This is described in Fig. 3. 

This factorization allows the use 
of common parameters to describe the 
three processes. The Virtuality of the 
process is given by Q^, the mass of the 
virtual photon, in inclusive and exclusive 
light vector meson DIS. In exclusive DIS 
production of heavy vector mesons this 
mass equals that of the produced meson 

(J/ip). For diffractive dissociation into two jets the mass of the di-jets is playing 
this role. From simple kinematics and assumimg that the masses of the jets are 
small compared with the mass of the di-jets, this can be shown to equal ^ 
kt is the transverse momentum of each jet. 
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Figure 2: measurement of the quark 
momentum wave function of the pion 
by measuring longitudinal momentum 
of the two jets from diffractive dissoci- 
ation. 
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Figure 3: Factorization of inclusive DIS, vector meson production and pion disso- 
ciation. 



2.3 Fermilab E791 Experiment 




Fermilab experiment E791 [6] recorded 
2 X 10^° events from interactions of 
a 500 GeV/c vr" beam with Car- 
bon (C) and Platinum (Pt) tar- 
gets. The trigger included a loose 
requirement on transverse energy 
deposited in the calorimeters. Pre- 
cision vertex and tracking informa- 
tion was provided by 23 silicon mi- 
crostrip detectors (6 upstream and , 
17 downstream of the targets), ten 
proportional wire-chamber planes, 
and 35 drift-chamber planes. Mo- 
mentum was measured using two 
dipole magnets. Two multicell, 
threshold Cerenkov counters were 
used for n, K, and p identifica- 
tion. Only 1/3 of the E791 data teraction of 500 GeV/c 
was used for the analysis presented Platinum targets. 
here. 
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Figure 4: qf distributions of di-jets from in- 

TT" with Carbon and 



The data was analysed by selecting events in which 90% of the beam momen- 
tum was carried by charged particles. This reduced the effects of the unobserved 
neutral particles. The selected events were subjected to the JADE jet-finding 
algorithm [7]. The algorithm uses a cut-off parameter {rricut) that describes the 
mass of individual jets. Its value was determined from Monte Carlo simulations to 
optimize di-jet detection. The di-jet invariant mass was calculated assuming that 
all the particles were pions. Only two-jet events were selected for further analysis. 
To insure clean selection of events, a minimum kt of 1.5 GeV/c was required. 

Diffractive di-jets were identified through the e^^''* dependence of their yield 
{qf is the square of the transverse momentum transferred to the nucleus and b is 
propotional to the r.m.s. radius squared of the nuclear target). Figure 4 shows 
the qf distributions of di-jet events from platinum (top) and carbon (bottom). 
The different slopes at low q^ in the coherent region reflect the different nuclear 
radii. Events in this region come from diffractive dissociation of the pion with the 
nucleus retaining its characteristic shape and remaining intact. 



2.4 Results 

The basic assumption that the momentum carried by the dissociating qq is trans- 
ferred to the di-jets, was examined by Monte Carlo (MC) simulation. Two samples 



of MC with 6 GeV/c^ mass di-jets were generated with different x dependences 
at the quark level. One sample simulated the Asymptotic wave function and 
the other the Chernyak Zhitnitski function. The two samples were allowed to 
hadronize through the LUND PYTHIA-JETSET [8] simulation and then passed 
through simulation of the experimental apparatus. 

In Fig. 5 the initial distri- 
butions at the quark level are 
compared with the final distribu- 
tions of the detected di-jets, includ- 
ing distortions in the hadroniza- 
tion process and influence due to 
experimental acceptance. As can 
be seen the qualitative features 
of the two distributions are re- 
tained. 



For events in the coherent re- 
gion described in the previous sec- 
tion the value of x was com- 
puted from the measured longitu- 
dinal momentum of each jet (eq. 
4). The resulting x distribution is 
shown in Fig. 6. Superimposed on 
the data is a fit to a combination of 
the two simulated final di-jet distri- 
butions of Fig. 5. 
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Figure 5: Monte Carlo simulations of the two 
wave functions at the quark level (left) and 
of the reconstructed distributions of di-jets 
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Figure 6: The preliminary x distribution of the diffractive di-jets fron the platinum 
target. The line is a fit to a wave function with about 90% asymptotic and 10% 
CZ. The CZ contribution is shown with a dashed line. 



The preliminary result under 
the analysis conditions as described 

is that the data can be fitted well with > 90% of the Asymptotic wave function 
and < 10% of the CZ wave function. The requirement that kt > 1.5 GeV/c can 
be translated to Q^ ~ 10 (GeV/c)^ (section 2.2). This shows that for these Q^ 
values the perturbative QCD approach that led to construction of the asymptotic 
wave function is reasonable. 

3 Pion DifFractive Dissociation and Color Trans- 
parency 

3.1 Theoretical Predictions 

DifFractive dissociation of pions to di-jets with kt > 1.5 GeV/c (Q"^ ~ 10 (GeV/c)^) 
can be associated with a \qq) system with size < r >\qq) ~ 1/Q < 0.1 fm. Such 
a small size configuration is expected to exhibit the effect of color transparency, 
namely a reduced interaction with the nuclear medium. The condition for ob- 
servation of the effect is that the system does not expand while traversing the 
nucleus. The distance that the system passes before it decays and expands (the 
"Coherence Length") is estimated by Ic ~ M^-m'^ where M is the di-jet mass. It 
is required that /^ > Ra- In E791 2piab = 1000 GeV/c. Using M ~ 6 GeV/c^ 
the result is Ic ~ 20 fm. which is larger than all nuclear radii. 

Bertch et al. [9] proposed that when high momentum pions hit a nuclear 
target the small \qq) component will be filtered by the nucleus and materialize as 
diffractive di-jets. Frankfurt et al. [10] proposed that diffractive dissociation of 
pions to high kt di-jets selects the small size \qq) Fock state with the size controlled 
by kt. They calculate the cross section for this diffractive dissociation of pions on 
nuclear targets with the striking result that for t = and kt > 1.5 GeV/c 



2.2 ...... 2.. „^' (^ = 0) (5) 



dxdMjd^qt dxdMjd'^qt 

Namely that the differential cross section at t = is proportional to A^ which 
is a very unusual dependence. The t-dependence for t 7^ is given by the nuclear 
form factor: |FA(t)|2 ~ e''*, b= ^^. 

3.2 Expected DifFractive Di-Jet A-dependence 

Realistically, it is impossible to measure the cross section at t = for two reasons: 
there is always a minimal t transferred to the target as a result of excitation of 
the pion to a heavy mass, and the statistical and systematic uncertainties near t 



= make it difficult to obtain precise results. 

The minimal t transfer can be estimated from simple kinematics. The result for 
a pion being excited from 771^^ to the di-jet mass Mjets is a longitudinal momentum 
transfer of: 

^^ = 2i. ^'^ 

with El the bombarding energy. The energy transferred to the nucleus is: 



For a numerical estimate we use: 



-t < 1/6 = 2 Carbon ^ 2.0 x lO'^ (GeV/c)^ (8) 

r^ fm 

Platinum ^ 4.3 x 10"^ (GeV/c)^ 

t X Ei/Ma ~ 0.83 (GeV/c)2 (carbon) and 0.011 (GeV/c)^ (platinum), ml = 0.02 
We can therefore neglect all the terms except for Mj^^^ and write 

Qz ^ ^^ (9) 

The minimal longitudinal momentum-square transferred to the nucleus as a 
consequence of exciting the pion to Mj^ts is q^ = —tmin- With M^efs ~ 6 GeV/c^ 
we get: 

-t™„ = ql - {^f - {^f = 0.0013 (GeV/cr. (10) 

While this value is small, it is not negligible compared with the t values for Pt. 
In order to avoid this problem and to increase the precision of the measurement 
we integrate the data in the diffractive region. This modifies the expected A 
dependence. We first write: 

^ = aj^A^e'' = ^^^26-^(^0+-^'+^?) (11) 

dt ^ ^ 

with qf the transferred transverse-momentum squared. From the numerical esti- 
mates we see that ^q can be neglected so that t = tmin—q^ sj[i(ida/dt -^ da/dqf. 
Integration of the diffiactive region (using R = RqA^^^) yields: 



{da/dq^)dq^ = a^A^e'*— e"''^* x (1/6) = a n A^ e''-- e~'"^^ x ——- (12) 

RqA^/-^ 



and the expected general A dependence is: ^^^ oc A"^^^ 

In the particular case of E791, using the numerical values of the nuclear radii 
of C and Pt: 

— ^ = =4.665= ; a; = 0.55. 13 

Re \2AA V 12/ ' ^ ' 



so that we will expect -^j^jj^ oc A^^^^. 



3.3 Analysis and Results 

To obtain the A dependence of the diffractive yield from the experimental results 
of Fig. 4 we must realize that these results are a combination of three elements: 
coherent nuclear diffraction, incoherent nuclear diffraction (interaction with indi- 
vidual nucleons) and smearing of both components by the experimental apparatus. 
In order to extract the coherent nuclear diffractive yield we simulate the coherent 
and incoherent qf distributions of the dissociation by assigning them the nuclear 
and the nucleon form factors, respectively. The simulated distributions are passed 
through the simulated experimental system so that the individual reconstructed 
distributions are obtained. At this stage a combination of the reconstructed distri- 
butions of the Pt nucleus and that of the nucleon is fitted to the data from the Pt 
target giving the ratio of coherent/incoherent yields and an overall normalization 
factor between the data and the simulation. The same is done for the C target. 
Results of these fits are shown in Fig. 7. 

The final results are obtained by integrating over the diffractive regions (of Pt 
and C) in the generated MC multiplied by the total normalization factor between 
MC and data. The analysis was done separately in two kf ranges since, as predicted 
by [10] there can be some dependence of a in a oc A" on kf. The PRELIMINARY 
result, compared with those expected according to [10], labelled a(CT) are: 



kt range a a(CT) 



1.5 < kt< 2.0 1.61 ± 0.08 1.45 

2.0 < kf 1.65 ± 0.09 1.60 



which are close to the Color Transparency expectations. 
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Figure 7: q^ distributions of di-jets for C and Pt targets. The lines are fits of the 
MC simulations to the data. 



4 Summary 



The following observations have been made: 



Diffractive dissociation of hadrons into di-jets can be used to study the 
hadron's internal quark structure. 

The momentum wave function of the \qq) in the pion is described well at 
Q^ ~ 10 (GeV/c)^ {kt > 1.5 GeV/c) by the Asymptotic wave function. 

Using data from Pt and C targets and assuming a oc A" the cross section 
for diffractive dissociation of pions to di-jets is found to be proportional to 
~ ^i.6i±o.o8 fo^ ^5 <: kt < 2.0 GeV/c and ~ ^i-esio.og ^^ 2.O GeV/c 

< kt . These results are consistent with predictions based on the color 
transparency effect. 
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